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We present the theory of a pair of atoms in a one-dimensional optical lattice interacting via a 
narrow Feshbach resonance. Using a two-channel description of the resonance, we derive analytic 
results for the scattering states inside the continuum band and the discrete bound states outside the 
band. We identify a Fano resonance profile, and the survival probability of a molecule when swept 
through the Bloch band of scattering states by varying an applied magnetic field. We discuss how 
these results may be used to investigate the importance of the structured nature of the continuum 
in experiments. 
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Ultracold atoms offer a compelling system for the study 
of fundamental physics, due to a large array of experi- 
mental tools allowing precision control of the state of the 
system. Two such tools are optical lattices [T] and Fes- 
hbach resonances 0. With an optical lattice the atoms 
are trapped in a periodic light potential, while a Fesh- 
bach resonance is used to tune the interactions between 
the atoms. Recently, these tools have made interesting 
two-body U U M H] and many-body [gj EU EE] 
states attainable. 

In this paper, we present an analytical description of 
the interplay between molecular binding and quantum 
motion in periodic potentials. We consider a pair of 
atoms in a one-dimensional optical lattice near a nar- 
row Feshbach resonance. Such one-dimensional geome- 
tries are realized in a three-dimensional optical lattice 
when the depth of the transverse cosine potential is much 
larger than the longitudinal lattice depth [3J U] . In a dis- 
crete lattice model [3] [5], we include on-site interactions, 
and consider the physics in the lowest Bloch band with 
nearest neighbor tunneling in the lattice. A two-channel 
model, in which a bound state in a closed channel is 
coupled to the continuum, describes our Feshbach res- 
onance [T5]. The formal description of our system 
bears a strong resemblance to Fano- Anderson theories of 
autoionizing states in atoms [T3] and localized electron 
states in solids |14j . 

The distinguishing feature of the lattice system is the 
structured nature of the continuum, which is divided into 
Bloch bands. As we are restricting our model to the low- 
est energy band, it means that the continuum has an 
upper edge. This has interesting and observable con- 
sequences on the interplay between bound and scatter- 
ing states. Two of these, the transmission Fano profile 
and the survival probability of a molecule when swept 
through the continuum by varying an applied magnetic 
field across the resonance, are explicitly derived, and we 
discuss how they may be probed in current experiments. 

Introducing center of mass, Z = (zi + Z2)/2, and rela- 
tive, z = zi — z 2 , coordinates the lattice allows a separable 



ansatz for the relative and center of mass motion: 



ip(zi, z 2 ) 



AKZ 



4>k(z). 
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The first Brillouin zone for the center of mass momentum, 
K, runs from —ir/a to n/a, where a is the lattice spac- 
ing. The Hamiltonian for the relative motion in the open 
channel, H°p = [-2JAf + E K + U5 Z . Q ] = H + U5 zfi , 
contains the discrete Laplacian 



K 



E_ 

Tj 



[ip K (z+a)+ip K (z-a)-2i) K (z)], (2) 



and the center of mass energy, Ek = — 4Jcos (Ka/2), in 
addition to the interaction between the atoms in the open 
channel. The on-site interaction in the entrance channel 
is U. We take the zero of energy to be the middle of the 
band, which spans the energy interval [— \Ek |, l-^ifl]- 

The interacting Green's function is found by explic- 
itly solving the Dyson equation, G K (E) — G° K (E) + 
G K (E)UG K (E). In coordinate space the result has the 
simple form 



G u K (E,z) = {z\G u K {E)\Q) 



1-UG° K (E,0) 



(3) 



where G° K (E,z) = (z\G° K (E)\0) is the Green's function 
for the relative motion of two non-interacting atoms on 
the lattice, which for energies inside the band (\E\ < 
\Er\) is propagating 



G° K (E,z) 



i exp(ip|z|) 



(4) 



with pa = cos 1 (E/Ek), while it falls of exponentially 
outside the band (|J5| > \E K \): 

G° K (E,z) = sign(£)^|t^ [-sign^)]* A\ (5) 



K 



Here na = cosh 1 \E/Ek\- For energies above the band 
the sign of the G° K (E, z) alternates between lattice sites. 
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The pole of G u K {E,z) at 1 = UG° K (E°,0) gives the en- 
ergy of the bound pair state; E% = 
For U > the pole constitutes the repulsively bound pair 
studied recently [3J- 

In the following we derive the resonant scattering prop- 
erties of the atoms and characterize the decay of molec- 
ular states in a two-channel model. We suggest avenues 
for the experimental investigation of these effects. 

Our aim is to solve the discrete coupled channel equa- 
tions for the relative motion on the lattice: 



EW*{z), (6a) 



H cl ^(z) + WS z ,o^(z) = E^(z), (6b) 

with a delta function coupling between the channels, 
(z,op\W\z', cl} = W5 Z , 5 Z , Z ,. 

From the spectrum of H cl we only retain the closed 
channel bound state, \4> rC s)i which causes the resonance. 
It is detuned by an energy E rcs (B,K) from the open 
channel band center and unit normalized, (4> rcs \(t>rcs) = 1- 
Hence we treat the closed channel Green's function, 
Gk(E) — [E — £Z" cl ] _1 , within a single-pole approxima- 
tion 



Gi(E) 



|</>rcs)(0r 



E-E Ies {B,K) 



(7) 



The bare resonance energy varies linearly with the ap- 
plied magnetic field, the slope being the difference be- 
tween the closed and open channel magnetic moments, 
A/i. To simplify the notation we introduce an effective 
coupling (z, op|Vy|(/> ros ) = WS Z: o. Both W and U may be 
obtained from free space scattering parameters by inte- 
grating over 3D tight-binding orbitals, and can be tuned 
by adjusting the transverse lattice. 

A numerical analysis based on perturbation the- 
ory shows that our restriction to the lowest band 
is valid provided the dimcnsionlcss coupling strength 
(amu/m)(ab g /100a )(A/z//iB)(A_B/G) is small com- 
pared with 0.01(^0 + 3)(lg" L )~ 6 , where m is the atomic 
mass, abg is the background scattering length, and AB 
is the magnetic field width of the resonance. The depths 
of the cosine optical lattice in the transverse, Vq , and 
longitudinal, Vq , directions are both measured in recoil 
units. The scaling factors are the atomic mass unit, amu, 
the Bohr radius, a , the Bohr magneton, ^b, and Gauss, 
G. If this criterion is not met, coupling to higher bands 
may not be neglected [XSl EE] ■ 

The scattering amplitude, f(E,K), is deduced from 
the open channel component of the asymptotic solu- 
tions to Eqs. ( |6a|6b |. It can be written as f(E,K) = 
fbg(E, K) + f Ies (E, K), where the background contribu- 
tion, f bs (E, K) = UG° K (E, 0)/[l - UG° K (E, 0)] is associ- 
ated with scattering in the absence of the coupling and 



has been described in [3J, while the resonant part is 
W 2 G° K (E,0)/[1-UG° K (E,0)} 2 



fres(E,K) 



E-E Yes (B,K)-GU(E,0)W 2 



(8) 



The scattering states cover the same range of energies, 
£if(fc) = Ek cos(fca), as the non-interacting pairs, and 
outside the band the poles of (l8|), E^, determined by 



El 



E 2 



Usiga(E b ) 



W 2 sign(E b ) 
Eb - E rcs (B, K) 



(9) 



correspond to discrete bound states of the coupled sys- 
tem. The bound state energies may be identified as a 
subset of the solutions to a quartic equation. The appar- 
ent pole at the uncoupled bound state energy, E®, has 
no weight, since the residues from the background and 
resonant parts of f(E, K) cancel. 

The size of the bare resonance state is much smaller 
than the lattice spacing, and thus confined to a single 
lattice site, the bare molecules effectively experience a 
lattice twice as deep as the one in which the atoms move. 
Hence the molecular tunneling rate is greatly suppressed 
with respect to its atomic counterpart, and we neglect in 
the following the modulation of the resonance energy by 
the center of mass motion, E Tes (B, K) « E ICS (B). 
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FIG. 1: (Color online) Scattering and bound states in an 
extended zone scheme for E lea = 2 J (a). The colormap inside 
the band illustrates the reflection coefficient, R, ranging from 
(white) to 1 (blue). The solid lines are the bound state 
energies of the coupled system, while the dashed lines indicate 
the diabatic, uncoupled bound states, with the horizontal line 
giving -Eres, and the other corresponding to E°. The on-site 
interaction and the coupling strength are U — 3J and W = 
1J, respectively, (b): closed channel fraction of the bound 



states, Eq. ( 11 1, (c): resonance profiles (solid lines) for E Tes = 
2 J at Ka/ir = 0, 0.3, 0.5, and 0.6 (from below). The dashed 
lines are the corresponding curves for the uncoupled case, 
where there is no resonance. 
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Figures [lja) and [2ja) show the variation of the scat- 
tering properties and the energies of the discrete bound 
states with the center of mass momentum and the bare 
resonance energy, respectively. The scattering is repre- 
sented through the reflection coefficient, R= |/| 2 , which 
ranges from (no scattering) to 1 (unitarity). The uni- 
tary scattering at the band edges is evidence of the van 
Hove singularities in the density of states, lmG° K {E, 0). 
In the uncoupled case, where W = 0, two diabatic bound 
states exist, as indicated by the dashed lines in Fig- 
ures [lj a) and[2]ja). The bare resonance state (horizontal 
line in Fig. [IJa), sloped line in Fig. [2ja)) gives rise to the 
clearly visible scattering resonance where it cuts through 
the band. The bound state pole, E®, of the open channel 
Green's function ^ is represented by the sinusoidal and 
horizontal dashed lines in Figures [lj a) and[2^a), respec- 
tively. For a finite coupling the two bound states repel 
each other, leading to the avoided crossing in Fig. [2j 

The coupling of the bare resonance state to the con- 
tinuum gives it a finite width, whenever it is embedded 
in the band. However, the resonance only exists for a 
range of center of mass motion Bloch states for a given 
magnetic field, as shown in Fig. [IT a) 7 . Outside this 
range (and between higher bands [17] ) there is instead a 
true bound state of the system. 





\ 













-5 5 10 15 

E les [J] 



-4-2 2 4 

E [J] 



FIG. 2: (Color online) (a): as Fig. [TJa), but plotted as a 
function of -E r es for fixed Ka/n = 0. The horizontal dashed 
line indicate the energy of the bound pair for W = 0, and 
the slanted dashed line gives the bare resonance energy, -E rC s- 

(b) : closed channel fraction of the bound states, Eq. flll) . 

(c) : resonance profiles for Ka/n = at E Tes /J = —2, -7, 5, 3 
(from the top at E/J — 0). 

The analytic shape of the resonance is described by a 
Fano profile in the transmission probability 



T =1 — R=T\ 



(10) 



arising due to interference between open and closed chan- 



nel scattering. The reduced variables e = 2(E — E les — 
A)/(hT) and q = 2A /HT contain shift and width func- 
tions related to the the real and imaginary parts of 
the molecular self energy: A(E,K) - iKT(E,K)/2 = 
{& CS \WG U K {E)W\^ CS ) = -W 2 {U+i^El-Ei)/{E 2 K + 
U 2 — E 2 ), while the background transmission coefficient 
is Tb g = 1/(1 + q 2 ). Note, that T assumes a simple 
Breit-Wigner profile, when the background interaction is 
ignored (U = 0, T bg = 1) [7J. 

We are aware that the lowest band approximation un- 
der which our calculations are valid is rather restrictive. 
However, to lowest order in the atom-molecule coupling 
the inclusion of higher bands only gives rise to real cor- 
rections to the molecular self-energy, modifying the shift, 
A, of the resonance energy, while leaving the lifetime un- 
changed. Furthermore, the shift will be only weakly en- 
ergy dependent, because the considered energies are well 
separated from the higher bands. The Fano resonance 
profile in the transmission is thus expected to be essen- 
tially unchanged if higher bands are included. 



The Fano profile (10 1 describes the scattering states 
of the system, but it also parametrizes the probability 
of transitions from a different initial state to the energy 
eigenstates inside the band. In the same way as Fano's 
theory was originally used to describe photoionization 
via autoionizing states [T3], Eq. (10) and Figures [TJc) 
and[2jc) describe the probability of RF- or photodissocia- 
tive transitions from deeper bound states in the molecule. 

We now turn to the bound states of the coupled sys- 
tem. Their wavefunctions have open and closed channel 
parts given by ^° b p (z) = yffiWG%(Et, z) and ^f(z) = 
\fZb5zfi, respectively. From normalization it follows that 
the closed channel fraction of the bound state at energy 
E b is 



Z h = 



W 2 \E b \ 

[l-UG° K (E b ,0)}^(E 2 -E 2 K )^ 



, (11) 



which is plotted in Figures [TJb) and[2^b). For the Fesh- 
bach molecule Z b vanishes exactly when the state enters 
the continuum, becoming the scattering resonance. 

In momentum space the component of the coupled 
channels bound state in the open channel is 



TO = 



2tt 1 



G° K (Eb,k) 
UG° K (E b ,z 



0)' 



(12) 



while (j>f{k) = ^aZ b /2ir for all k. The non- interacting 
quasi-momentum Green's function is Q° K {E,k) = [E — 
E K cos(ka) + iri\^ 1 . For \K\ < n/a and for E b > (E b < 
0) the bound state function 4>°^{k) is peaked at the edges 
(center) of the Brillouin zone. This reflects that <j>°^(k) 
is predominantly comprised of the fc-states with energies 
in proximity of E b . For E b > (E b < 0) this corresponds 
to the states near the top (bottom) of the band [S] . Thus 
the location of the peaks makes it possible to distinguish 
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between bound states lying above or below the continuum 
in time of flight [3] . The height of the peaks are reduced 
as K approaches ±7r/a, where all quasi-momentum states 
are degenerate in our model. 

At a given energy inside the band the Fcshbach 
molecules are unstable. Their dissociative decay width 
is HT(E, K) given above. For a sweep of the magnetic 
field across resonance, starting with the resonance en- 
ergy tuned below the band, the molecular population 
follows a simple rate equation, dm/dt = —T(E(t))m(t), 
provided the sweep is fast enough that recombina- 
tion may be neglected. For a linear sweep, where 
dE/dt — A/iB, the remaining molecule fraction can be 
found analytically as a function of the energy inside 
the band: m(E) = exp{-2a[(sin" 1 E/\E K \ + tt/2) - 
C^an- 1 CE/\E K y\ - (E/\E K \) 2 + tt/2)]}, where a = 
W 2 /AfiBh, (5 = U/E K , and C = y^JJl+Wj. A 
unique feature of the lattice system, is that the contin- 
uum has an upper edge, and for a sweep ending above 
the band the final molecule fraction has the simple form 



X = m(\E K \) = cxp 



-2ira 1 - 



P 2 



1 + /9 2 



(13) 



depending only on the two constants a and (3. This ex- 
pression is plotted in the left hand panel of Fig. [3] for 
U = 2J. Again, this result is expected to be modified 
only slightly by the inclusion of higher bands. 

In the case of vanishing background interaction, U = 0, 
the survival probability assumes a Landau-Zener form, in 
which a plays the role of the adiabaticity parameter. We 
also note that in the limit U — » ±oo all molecules survive 
the passage through the band. The limit a — > 0, corre- 
sponding to sweeps instantaneous on the time scale set 
by the coupling, also leaves the molecules intact. Finally, 
X — > 1 at K = n/a due to the vanishing width of the band 
at the center of mass zone boundary. This is an artifact 
of omitting tunneling beyond nearest neighbor. As more 
distant tunneling amplitudes are incorporated the band 
widens at the zone boundary and an extra van Hove sin- 
gularity develops in the density of states [5] . Lastly, we 
note that x is independent of the sign of U. 

The average energy of the dissociated atoms is plotted 
in the right hand panel in Fig. [3} As a smaller fraction of 
the molecules survive the passage through the band the 
dissociation energy per atom approaches half the lower 
band edge, — \E K \/2 — — 2J cos(Ka/2). Molecule disso- 
ciation energies have been measured in free space [18|H9]. 
and the energy distribution after the sweep through the 
energy band may be identified by studying the ballistic 
expansion of the gas after the atoms are released from 
the lattice on a time scale, admitting adiabatic mapping 
of the quasi-momentum states to plane wave momentum 
states pQ. 

N. N. acknowledges financial support by the Danish 
Natural Science Research Council. 
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FIG. 3: Remaining molecule fraction after a linear sweep 
through the band (left) and average dissociation energy per 
atom (right) for a=0, 0.1, 0.5, 1, 2, 5, 10 (from the top). 
U — 2J. For fast ramps all molecules survive, but as the 
ramp speed is decreased more molecules decay and the disso- 
ciation energy approaches half the lower band edge (dashed 
line). 
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